Generation and Measurement of Non Equilibrium Spin Currents in Two Terminal 

Systems 
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Generation and measurement of non-equilibrium spin current in two probe configuration is dis- 
cussed, ft argued and shown that spin current can be generated in two terminal non-magnetic 
system. Further it is shown that these spin currents can be measured via conductance in two probe 
configuration when the detector probe is ferromagnetic. 
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The generation, manipulation and detection of spin 
currents, a flow of angular momentum, is a major goal 
of spin-based electronics. Thus establishing methods for 
efficient generation and detection of spin currents is a key 
for further advancement of spintronics. Amongst the dif- 
ferent approaches for spin current generation and manip- 
ulation, spin-orbit(SO) coupling is attracting consider- 
able interest. It has been predicted that in a SO coupled 
system a non zero spin current flows in a direction per- 
pendicular to the applied electric field [l[ . This gives rise 
to related spin-Hall effect (SHE) causing spin accumula- 
tion at the edges of sample 0, H[ . The spin accumulation 
caused by SHE has been detected by optical techniques 
in a semiconductor channel Q and in two-dimensional 
electron gas (2DEG) in semiconductor heterostructures 
0, [fjj . Recently electronic detection of spin accumulation 
in Hall cross geometry which is a four probe system, was 
reported in Ref. @, Q ■ Till date most of the theoretical 
and experimental studies are concerned with the obser- 
vation of SHE which is a spin accumulation caused by the 
flow of spin current in multi-probe systems. To the best 
of our knowledge the possibility of spin current genera- 
tion and detection in two probe systems via conductance 
measurement has not been reported yet. In this article 
we address this possibility. We consider two configura- 
tion shown in Fig. 1(a) and 1(b) where a 2DEG with SO 
coupling is sandwiched between two ideal leads without 
SO coupling. In Fig. 1(a) both leads are non-magnetic 
while in Fig. 1(b), the detector lead is ferromagnetic (FM) 
and injector lead is non-magnetic. 
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FIG. 1: Setup for spin current generation and detection (a) 
current leads are non-magnetic and injected currents is unpo- 
larized but outgoing current is polarized due to SO region (b) 
Detector lead is magnetic; and the conductance depends on 
absolute director of magnetization {6,<f>) because of the spin 
current generated by SO region 



An intuitive understanding of the non-equilibrium spin 
current generation can be gained by considering the spe- 
cific case of Rashba SO coupling which in 2DEG has the 
form, Hr = \{k x a v — k y o~ x ), where k x and k y are compo- 
nents of momentum vector in the plane of 2DEG and o x 
and <jy are Pauli matrices. Consider the situation when 
the incident current flows along positive x direction (Fig. 
1(a)). In this particular case the Rashba interaction is 
effectively, X(k x a y ) because k y =0. Therfore the outgo- 
ing current gets polarized along the y axis. In general 
SO interaction polarizes scattered beam. This general 
result is valid for arbitrary form of SO interaction viz. 
Rashbaj3jDresselhaus[lO| or impurity induced SO inter- 
action [J [Tl|, EH • However the direction and magnitude 
of polarization for the scattered beam depends on the 
exact form of SO interaction [13(. Therefore the outgo- 
ing charge current in non-magnetic lead two is polarized 
even though the incident current is unpolarized. Hence 
it is a true flow of non- equilibrium angular momentum, 
from SO coupled region to the detector lead. 

Let us now focus on the conceptual problem of non- 
equilibrium spin current measurement via conductance 
in two probe system with ferromagnetic (FM) detector 
shown in Fig. 1(b). According to Onsager's relation the 
conductance of a two terminal system with one contact 
being ferromagnetic is invariant upon magnetization re- 
versal, ie.,G(m)=G(— m), where rh is unit vector along 
the magnetization direction. In absence of SO interac- 
tion, system shown in Fig. 1(b) is rotationally invariant, 
therefore conductance does not depend on the direction 
of magnetization and Onsager symmetry is satisfied triv- 
ially. In presence of SO interaction since rotation symme- 
try is broken therefore conductance depends on the ab- 
solute direction of magnetization for configuration shown 
in Fig. 1(b). This is because the detector lead is ferro- 
magnetic (Fig. 1(b)), therefore electrons polarized in the 
direction of magnetization will transmit easily compared 
to the electrons polarized along another direction. Hence 
the conductance of NM/SO/FM system (Fig. 1(b)) will 
depend on the absolute direction of magnetization 9 and 
<f> and will be symmetric in 9 and <f> consistent with On- 
sager symmetry.In other words, G(rhi)^ G(rfi2), where 
rhi and rn.2 are two arbitrary unit vector along which 
magnetization is pointing, with the condition that rh\ 



7^ rhi- We show that these conductance variation as 
a function of absolute direction of magnetization is in 
phase with the polarization of outgoing current, there- 
fore, it measures the non-equilibrium spin current. More- 
over we show that to detect spin currents magnetization 
reversal is not necessary, rather tilting the magnetization 
away from the original direction is enough. Magnetiza- 
tion tilting can be achieved using external magnetic fields 
as is done usually [HJ|. The spin currents generated by 
SO interaction can also tilt the magnetization because it 
will give rise to a torque on ferromagnetic detector. The 
torque arises because the spin current generated by SO 
interaction in general will have components transverse 
to the magnetization direction of Ferromagnetic contact. 
The transverse component of spin current is transferred 
to the FM as torque [15J. If the torque is large enough 
it can tilt the magnetization. For nm scale thick fer- 
romagnetic film it is possible increase the torque by in- 
creasing the current density 16]. In case torque is not 
large enough to tilt magnetization one can always use an 
external magnetic filed to rotate the magnetization. Re- 
cently Adagideli, Bauer and Halprerin in Ref. [l7| have 
also discussed the possibility of detecting torque related 
signal via small angle magnetization dynamics in three 
terminal systems. We stress that our proposal works in 
two terminal systems and relies on conductance measure- 
ments which are easier as well have potential to be used 
in already developed semiconductor technology. 

To study the above discussed effect quantitatively we 
make use of spin density matrix scattering theory de- 
veloped in Ref. [l5[ which treats the non-magnetic and 
magnetic system at the same footing. An incident cur- 
rent in channel \n,a) = \n) ® \a)(n channel index due 
to transverse confinement along y direction and a is spin 
index) leads to the following spin density matrix for out- 
going current (l5| . 
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where N is normalization factor ensuring that Tr(p) = 1 
and and are spin resolved transmission coef- 

ficient. The incident current is unpolarized (completely 
mixed state) therefore, the spin density matrix for the in- 
cident current is pi n = n a \a) + n_ Q | — a) with n Q =n_ Q . 
Here n a and n_ a is the number of up electrons and down 
electrons respectively. The fact that incident current is 
unpolarized leads to the following spin density matrix for 
outgoing current, 
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From eq[2] and eqfT] it is clear that spin density ma- 
trix for outgoing current is completely determined if spin 
resolved transmission coefficients are known. Since the 
polarization direction of outgoing current in general will 
lie along any arbitrary direction, therefore, the transmis- 
sion coefficient should be calculated in eigen-basis of <x • u 




FIG. 2: (Color online) Spin polarization of outgoing cur- 
rent versus quantization axis (0=90 and (j> is varied) for 
non-magnetic systems corresponding to fig. 1(a). The di- 
mensionless SO interaction strength a=0.1 and Ef =l-0t and 
kfl — 100a(l is mean free path) for diffusive case. Straight 
line with filled circle corresponds to the magnitude of polar- 
ization. 



where u = (sin cos 0, sin#sin</>,cos#) is the quantiza- 
tion axis. The polarization along the quantization axis 
u is given by the trace of matrix product of <x • u and 
Pf given in cq[2J However to determine polarization vec- 
tor we need to calculate polarization components along 
two more direction which are perpendicular to u and are 
linearly independent. Therefore, for simplicity and clar- 
ity we choose two specific cases of quantization axis u, 
namely (a) We fix 9 — it/2 and vary <f>. In this case 
three linearly independent axis are h\ = (coscf), sin<^>, 0), 
h 2 = (— sin 4>, cos </>, 0) and h 3 = (0, 0, 1) = z. (b) We fix 
<f> = and vary 9. In this case three linearly independent 
axis are hi = (sin 6*, 0, cos 9), h 2 — (— cos 9, 0, sin#) and 
hi = (0, 1, 0) = y. The polarization along a specific axis 
hi is calculated using, 
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1,2,3 
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where pf is density matrix given by eqU] 

To obtain results we perform numerical simulation for 
a two dimensional electron gas (2DEG) using spin de- 
pendent recursive Green-function technique. The details 
of this can be found in Ref. [ll|, [l2j]. We model the 
conductor on a square tight binding lattice of size 50 x 50 
with lattice spacing a and hopping parameter t. We fix 
the Fermi energy Ep = 1.04 and take t=l as the unit of 
energy. The tight binding dimensionless SO parameter 
is a = \/(ta). For disorder we take Anderson model, 
with on-site energies distributed randomly between [- 
U/2, U/2], where U is the width of distribution. The 
results presented are for Rashba Spin-Orbit coupling 0] . 
However the conclusion remain valid for other type of SO 
couplings as well. 

The three linearly independent polarization compo- 
nents of outgoing current are shown in Fig. 2 for non- 
magnetic configuration shown in Fig. 1(a). We have fixed 
9 = 7r/2 and rotated the quantization axis in xy plane 
which correspond to the case (a) discussed above. For 



3 




,i — , — i — , — i — , — i — , — i i — , — i — , — i — , — i — , — 1„ 2 

90 180 270 360 90 180 270 360 

9 or (p(degree) 9 or (p (degree) 



FIG. 3: (Color online) Conductance and Spin polarization of 
outgoing current as a function quantization axis for ballistic 
system with ferromagnetic detector shown in Fig. 1(b). Other 
parameters are same as in fig. 2. 

ballistic system (Fig. 2, left panel) outgoing currents is 
perfectly polarized along y axis (curve for Pi has maxi- 
mum value at 0=0, and 2n while P\ has maximum and 
minimum values at cf> — tt/2 and 3ir/2 respectively.) This 
is because current is flowing along positive x axis the ef- 
fective Rashba interaction becomes X(k x a y ) as mentioned 
in the introduction. However in the presence of disorder 
due to scattering the current direction in 2DEG changes 
randomly therefore the Rashba coupling retains its two 
dimensional nature. Hence all three components of po- 
larization are non zero for diffusive system (Fig. 2, right 
panel). Moreover the magnitude for each component as 
well the net polarization (shown by straight line with 
filled circles in Fig. 2) is large in diffusive regime com- 
pared to the ballistic case. This is because in presence 
of disorder electron mobility decrease which in turn in- 
crease the effect of Rashba coupling. Therefore moderate 
disorder will be helpful in generating spin current which 
is promising from the experimental point since it will not 
require very clean or ballistic sample. 

From Fig. 2 it is clear that the outgoing current is spin 
polarized though the leads are non-magnetic (Fig. 1(a)). 
Since the leads are non-magnetic therefore the transmis- 
sion coefficient will be independent of spin polarization 
of outgoing current and as a consequence conductance 
will also be independent of quantization axis. Hence with 
non-magnetic leads one cannot measure the spin current. 
However, if the detector lead is magnetic (configuration 
shown in Fig. 1(b)) then transmission coefficient will be- 
come spin dependent which will cause conductance to 
depend on quantization axis as discussed in paragraph 
3. This is confirmed in Fig. 3, where we have plotted the 
conductance(left panel) and polarization (right panel) for 
the configuration corresponding to Fig. 1(b). Ferromag- 
net detector is modeled as exchange split with exchange 
splitting (A) given as A/Ep=0-5 The quantiza- 

tion axis is given by the direction of magnetization. We 
have plotted results for two different quantization axes 
discussed in paragraph 4 above, namely, (1) 4> =0 and 9 
is varied, rotation in xz plane (2) 0=90 and <fi have been 
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FIG. 4: (Color online) Conductance and Spin polarization of 
outgoing current a function of quantization axis for diffusive 
system with ferromagnetic detector shown in Fig. 1(b). The 
results are for diffusive case with kfl = 100a(Z is mean free 
path). Other parameters are same as in fig. 3 



varied, rotation xy plane. We notice that conductance 
and polarization for zero SO coupling is independent of 
quantization axis shown by straight line with filled circles 
in Fig. 3. A finite polarization in absence of SO coupling 
is due to the presence of ferromagnetic lead. In presence 
of SO interaction polarization increases compared to the 
polarization for zero SO coupling which means spin cur- 
rents is being generated by SO coupling. Moreover po- 
larization and conductance show in phase variation as a 
function of 9 or 0. Therefore conductance measures the 
spin current as discussed before. We notice that, the po- 
larization curve in Fig. 3 for 9=90 have different values 
for 0=90 and 0=270 which corresponds to y and -y axis 
respectively. This is because in ballistic system outgoing 
current is polarized along positive y axis as discussed ear- 
lier, therefore, the net polarization decreases (increases) 
when ferromagnet points anti-parallel (parallel) to y axis. 
The asymmetry in polarization can be reduced by reduc- 
ing the magnetization of FM (decreasing A/Ep) and by 
using a weakly coupled detector lead(non-invasive). This 
would allow one to measure spin currents just due to SO 
coupling. These finding are robust and survives in pres- 
ence of disorder as shown in Fig. 4. In presence of disor- 
der due to scattering the outgoing current gets polarized 
along some arbitrary direction therefore the asymmetry 
is visible throughout the 9 or range. However the po- 
larization variation still remain in phase with the conduc- 
tance. We stress that to measure spin current one does 
not need to reverse the magnetization because the con- 
ductance changes as soon as the magnetization is tilted 
away from the original direction as seen in Figs. 3 and 4. 

The magnetization can be tilted due to the torque 
which arises because of spin currents as discussed in 
paragraph three. Dimensionless torque is defined as[l5| 
r = ^§777^7' w here M is magnetization vector of FM, 
P is polarization vector of outgoing current and / is cur- 
rent. We have plotted the magnitude of torque in Fig. 5 
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FIG. 5: (Color online) Torque acting on detector ferromagnet 
as a function of quantization axis corresponding to configura- 
tion shown in Fig. 1(b). 

for ballistic (left panel) and diffusive system(right panel) 
respectively. For ballistic system the torque goes to zero 
when magnetization is either parallel or antiparallel to y 
axis (see Fig. 5 left panel, 9=90 and cf> variable) because 
outgoing current is polarized along y axis as discussed 
above and as a consequence torque vanishes. For all other 
magnetization direction and for diffusive system (Fig. 5, 
right panel), the torque is always non zero. If the cur- 
rent density is large enough and Ferromagnetic film is of 
nm scale thick then, magnetization can be tilted just by 



increasing current and as a consequence one would see a 
change in conductance which measures the spin current 
as mentioned earlier. 

Finally let us discuss the possibility of observing the 
spin current related effects presented above. From Figs. 3 
and 4 it is seen that conductance changes by 2-3% when 
magnetization direction is rotated. Magneto-resistance 
experiments can easily detect changes of this order [lj| • 
The 2-3% conductance change was obtained for dimen- 
sionless Rashba SO parameter a=0.1. This corresponds 
to the Rashba spin splitting energy of the order of lOmeV. 
The Rashba SO splitting of 30meV has been reported in 
Ref. [HI in II-IV HgTe system. Therefore, the dependence 
of conductance on absolute direction of magnetization 
can be easily detected which measures the spin currents. 
In addition, tunability of Rashba coupling via external 
gate voltage [l8j provides further control over the spin 
current generation, torque and associated conductance 
change. 

In conclusion, we have presented a general method to 
generate and measure non-equilibrium spin currents in 
two probe configuration via conductance. We have clar- 
ified that magnetization reversal is not required to mea- 
sure spin currents and shown that it is consistent with 
Onsager symmetry. Further we have show that torque 
arising due to spin current can be used to tilt the mag- 
netization which will help in measuring spin currents. 
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